UNCLASSIFIED 


_ AD  NUMBER _ 

AD838850 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
and  their  contractors;  Critical  Technology;  JUN 
1968.  Other  requests  shall  be  referred  to  Naval 
Postgraduate  School ,  Code  023 ,  Monterey,  CA 
93940.  This  document  contains  export-controlled 
technical  data . 


_ AUTHORITY 

usnps  ltr,  21  jan  1972 


THIS  PAGE  IS  UNCLASSIFIED 


UNITED  STATES 

NAVAL  POSTGRADUATE  SCHOOL 


mm 


THESIS 


A  METHODOLOGY  FOR 

INVESTIGATING  THE  FACTORS 

AFFECTING  Till 

OUTCOME  OF  A  SUBMARINE 

VERSUS  S 

UBMARINF  ENCOUNTER 

k  •  r- 

by  \  .  '  ■ ' 

O 

v . 

Ralph 

Edward  Tuggle 

\  >  '  .  f, 

\  ‘  c^°  ; 

_ u-U  ■■ 

June  I96S 


This  document  is  subject  to  special  export  con¬ 
trols  and  each  transmittal  to  foreign  government 
or  foreign  nationals  may  be  made  only  with  prior 
approval  of  the  U.  S.  Naval  Postgraduate  School. 

[(Code-  023,  Monterey  Calif. 


A  METHODOLOGY  FOR  INVESTIGATING  THE  FACTORS 
AFFECTING  THE  OUTCOME  OF  A  SUBMARINE 
VERSUS  SUBMARINE  ENCOUNTER 


by 


Ralph  Edward  Tuggle 

Lieutenant  Commander,  United  States  Navy 
B.S.,  Naval  Academy,  1959 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  1968 


Signature  of  Author  (/ ' 


Approved  by 


The-yis  A3vT 


sor 


f  Academic  Dean 


ABSTRACT 


A  quantitative  analysis  of  encounters  between  two 
submarines  was  made  in  order  to  investigate  the  effects 
different  factors  have  on  the  outcome  of  such  encounters. 
A  model  was  constructed  with  the  use  of  a  logic  tree  and 
a  computer  simulation  of  a  submarine  versus  submarine  en¬ 
counter.  The  computer  simulation  was  developed  from  the 
logic  tree.  The  outcomes  of  the  simulation  were  analysed 
using  contingency  table  tests  of  independence,  the  theory 
of  games  of  timing  and  a  linear  statistical  model.  The 
contingency  table  tests  and  the  theory  of  games  of  timing 
demonstrated  the  relationship  of  range  and  elapsed  time 
to  the  possible  outcomes.  The  linear  statistical  model 
was  used  to  obtain  estimates  of  the  effects  that  various 
own  ship  and  target  capabilities  have  on  the  outcome  of 
the  encounters. 
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CHAPTER  I 


INTRODUCTION 


Background 

Improvement  in  combat  force  capabilities  and 
effectiveness  has  been  and  continues  to  be  one  of  the 
major  objectives  of  the  military  planner.  In  fact,  since 
the  advent  of  sophisticated  and  costly  weapons  systems 
capable  of  practically  unlimited  destruction,  this  objec¬ 
tive  has  assumed  even  more  importance.  Trial  and  error 
methods  of  determining  needed  improvements  are  costly  and 
inefficient.  Proposed  changes  could  yield  little  improve¬ 
ment  in  capability  at  a  high  cost  whereas,  at  the  same  or 
perhaps  a  lesser  cost,  much  larger  gains  could  be  obtained 
in  other  areas.  It  has  long  been  realized  that  if  encoun¬ 
ters  between  opposing  forces  could  be  formulated  into  a 
mathematical  model  then  an  analysis  of  the  parameters  con¬ 
stituting  the  model  could  point  out  the  areas  in  which  im¬ 
provement  is  most  needed  and  would  be  most  meaningful. 

Previously  Developed  Models 

Several  types  of  models  have  been  developed  and 
utilized  in  the  past  and  are  still  evolving.  Some  of  the 
more  notable  of  these  are  the  Lanchester's  Equations,  game 
theory  (in  particular,  games  of  timing  or  stochastic  duels), 
war  gaming,  Monte  Carlo  simulations  and  other  manual  or 
computerized  simulations. 


TRW  Systems,  Incorporated,  of  Redondo  Beach,  Cali¬ 
fornia,  utilized  a  model  describing  an  encounter  between  a 
submarine  and  a  surface  ship  task  force  in  order  to  inves¬ 
tigate  certain  surface  ship  sonar  systems  and  ASW  weapon 
systems.  The  model  proposed  by  this  thesis  has  been  devel¬ 
oped  along  the  same  lines  as  the  TRW  simulation,  but  dif¬ 
fers  in  several  respects.  The  most  obvious  difference  is 
that  this  model  describes  an  encounter  between  two  sub¬ 
marines.  Whereas  the  real  elapsed  time  of  the  encounter 
was  important  in  the  TRW  investigation,  it  is  not  perti¬ 
nent  to  this  one.  TRW  was  specifically  interested  in  com¬ 
paring  various  sonar  and  weapon  systems  while  this  paper 
has  attempted  to  investigate  other  factors  as  well.  Fi¬ 
nally,  while  the  TRW  model  utilized  a  manual  evaluation 
of  the  logic  tree,  this  thesis  has  adapted  the  logic  tree 
evaluation  to  a  computer  simulation. 

The  Submarine  Versus  Submarine  Problem 

In  order  to  construct  a  valid  model,  it  is  neces¬ 
sary  to  have  an  understanding  of  the  different  factors 
which  compose  or  affect  the  progress  of  the  encounter. 

The  fire  control  problem  for  a  submarine  versus  submarine 
encounter  is  very  complex.  Achieving  an  accurate  solu¬ 
tion  is  much  more  difficult  than  in  a  context  where  the 
antagonists  can  "see"  each  other  either  electronically  or 
visually,  thus  making  relatively  accurate  range  informa¬ 
tion  an  integral  part  of  the  problem.  As  new  weapons 
that  are  capable  of  extended  ranges  and  new  sonars  having 


long  range  detection  and  classification  capabilities  are 
developed,  the  doctrine  and  equipment  utilized  in  obtain¬ 
ing  a  submarine  fire  control  solution  are  continually  being 
refined  and  modernized.  The  standard  method  of  obtaining 
a  solution  is  the  bearings-only  approach.  In  this  type  of 
an  approach,  bearing  rate  and  estimated  target  speed  are 
used  to  arrive  at  a  best  estimate  of  target  course,  speed 
and  range.  The  bearings  only  approach  has  the  character¬ 
istic  that  the  accuracy  of  the  solution  is  always  in  doubt 
although  it  tends  to  improve  with  an  increase  in  the  time 
that  the  contact  is  held. 

Obviously,  the  elapsed  time  between  initial  detec¬ 
tion  of  a  target  and  launching  a  weapon  will  affect  the 
outcome  of  the  encounter.  Current  doctrine  specifies  that 
in  order  for  a  submarine  to  most  effectively  use  its  weap¬ 
on  and  sonar  capabilities  and  at  the  same  time  maximize 
its  survivability,  it  should  shoot  at  as  long  a  range  as 
possible.  There  are,  of  course,  many  factors  other  than 
just  firing  doctrine  which  contribute  to  the  determination 
of  the  final  results  of  an  encounter.  The  sophistication 
and  reliability  of  the  weapon  system,  sonar  system  and 
fire  control  system  are  extremely  important  as  well  as  the 
state  of  training  and  experience  level  of  the  crew  in  uti¬ 
lizing  these  systems.  The  design  of  the  ship,  its  propul¬ 
sion  and  control  systems,  will  also  have  a  marked  impact 
on  the  results. 
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A  snap  shot,  as  the  term  is  used  in  this  paper,  is 
basically  a  defensive  tactic.  It  is  employed  when  a  ship 
has  been  attacked.  In  -such  a  situation,  a  homing  torpedo 
is  launched  in  the  direction  from  which  the  attack  was 
launched  or  is  expected. 

Organization  of  the  Remainder  of  the  Thesis 

The  objective  of  this  thesis  is  to  propose  a  meth¬ 
odology  for  the  investigation  of  the  effects  each  of  the 
various  factors  has  on  the  outcome  of  a  submarine  versus 
submarine  encounter.  The  fulfillment  of  this  objective  is 
accomplished  by  constructing  a  model  which  makes  use  of  a 
logic  tree  and  computer  simulation.  These  are  discussed 
in  detail  in  Chapter  II,  along  with  all  simplifying  as¬ 
sumptions.  The  evaluation  of  the  simulation's  results  are 
contained  in  Chapter  III.  Contingency  tables,  the  theory 
of  stochastic  duels  and  a  linear  statistical  model  were 
used  in  that  evaluation.  The  conclusions  reached  after 
the  analysis,  and  recommendations  for  further  study  are 
contained  in  Chapter  IV.  The  appendices  contain  the  var¬ 
ious  diagrams,  computer  programs  and  outputs  relating  to 
this  study. 
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CHAPTER  II 


APPROACH  TO  THE  PROBLEM 

The  initial  step  in  the  accomplishment  of  the  ob¬ 
jective  of  this  paper  was  to  build  a  valid  model  of  an 
encounter  between  two  submarines.  A  logic  tree  was  con¬ 
structed  to  graphically  represent  the  encounter.  From 
that  logic  tree,  a  computer  program  was  developed  to  sim¬ 
ulate  the  encounter.  The  computer  was  utilized  to  provide 
a  large  sample  number  of  encounters  in  order  that  a  more 
meaningful  analysis  of  the  factors  involved  could  be  made. 

The  Logic  Tree 

The  logic  tree  as  illustrated  in  Appendix  A  is  a 
connected  directed  graph  composed  of  probabilistic  nodes. 
The  nodes  represent  the  opportunity  for  the  occurrence  of 
the  pertinent  events  constituting  the  encounter.  Each 
event  has  a  positive  probability  of  occurring.  Each  node 
has  two  arcs  incident  from  it.  One  arc  represents  the 
actual  occurrence  of  the  event  while,  if  it  does  not  take 
place,  ''he  other  arc  is  followed.  For  this  model,  the 
occurrence  of  an  event  was  assumed  to  be  a  Bernoulli  tri¬ 
al.  The  path  the  encounter  follows  down  the  tree  is  de¬ 
termined  by  comparing  the  probability  associated  with  the 
first  node  with  a  random  number  drawn  from  a  uniform  (0,1) 
distribution.  If  the  random  number  is  less  than  or  equal 
to  the  nodal  probability  the  event  is  said  to  have  oc¬ 
curred  and  the  appropriate  arc  is  followed  to  the  next 
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node.  Likewise,  if  the  random  number  is  greater  than  the 
nodal  probability  the  event  did  not  occur  and  the  oppo¬ 
site  arc  is  followed.  This  procedure  is  repeated  until 
the  path  reaches  its  final  node  at  which  time  the  encoun¬ 
ter  is  terminated. 

To  aid  in  the  analysis  and  evaluation  of  the  out¬ 
comes,  the  tree  was  divided  into  two  cases.  Case  one 
represents  an  encounter  in  which  the  target  does  not  de¬ 
tect  the  presence  of  own  ship  until  after  own  ship  has 
launched  a  weapon.  Case  two  is  the  converse  of  case  one 
in  that  the  target  first  detects  own  ship  prior  to  own 
ship  launching  a  weapon. 


Definition  of  Events  and  Outcomes 

PCDET  -  Probability  target  detects  own  ship  prior 
to  own  ship  launching  a  weapon. 

PATT  -  Probability  own  ship  makes  a  successful 
attack  on  the  target. 

PCA  -  Probability  target  will  counterattack  prior 
to  own  ship's  attack. 

PCAT  -  Probability  the  target's  counterattack  is 
successful . 


PTAS 


PTSS 


PRAT 


-  Probability  target  attempts  a  snap  shot, 
given  the  target  has  been  fired  upon. 

-  Probability  the  target's  snap  shot  is 
succe  ssful . 

-  Probability  that  own  ship  makes  a  success¬ 
ful  snap  shot,  given  that  own  ship  has  been 
fired  upon. 

-  Probability  own  ship  makes  a  successful  re¬ 
attack  on  the  target. 


PF’VADl  -  Probability  target  evades  (breaks  off  the 
encounter)  given  he  has  been  attacked. 


PEVAD2  -  Probability  target  evades  given  he  detects 
own  ship  prior  to  being  attacked. 

DPAT  -  Probability  of  successful  attack  on  the  tar¬ 
get  given  the  target  detects  own  ship  prior 
to  being  attacked. 

UPl'S  -  Probability  target  makes  a  successful  snap¬ 
shot  given  the  target  detects  own  ship  prior 
to  being  attacked. 

Each  path  in  the  tree  terminates  with  one  of  four 
possible  outcomes  of  an  encounter.  These  outcomes  are: 

IW  -  Target  is  sunk  and  own  ship  survives  the  en¬ 
counter  (hereafter  referred  to  as  a  success¬ 
ful  attack). 

IL  -  Own  ship  is  sunk  and  the  target  survives  the 
encounter . 

IDW  -  Both  ships  survive  the  encounter. 

IDL  -  Both  ships  are  sunk. 

Computer  Simulation 

The  flow  chart  from  which  the  computer  program  was 
written  is  contained  in  Appendix  B.  If  was  written  in 
FORTRAN  IV  and  run  on  an  IBM  360  series  computer.  The 
computer  program  comprises  Appendix  C.  The  total  time  for 
compilation  and  execution  was  about  2.7  minutes.  During 
this  time,  a  total  of  twenty-seven  different  combinations 
of  own  ship  and  target  capabilities  were  involved  in  one 
thousand  encounters  for  each  of  three  regions,  yielding  a 
total  of  81,000  encounters. 

Parameters  and  Associated  Capabilities 

the  values  assigned  to  the  nodal  probabilities  re¬ 
flect  the  relative  capabilities  of  the  submarines  in- 
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volved.  By  simply  altering  the  parameter  values,  the  model 
can  be  made  to  describe  encounters  between  any  two  types  of 
submarines.  The  values  tor  PCDET  are  a  measure  of  the  tar¬ 
get's  sonar  system  capability  while  the  values  assigned  to 
PEVADl  and  PEVAD2  are  measures  of  the  target's  ship  system 
capability.  These  parameters  can  also  be  thought  of  as 
containing  some  yardstick  of  own  ship's  capabilities.  As 
an  example,  a  large  value  for  PCDET  might  indicate  a  noisy 
attacking  ship  instead  of  a  particularly  sophisticated  tar¬ 
get  sonar  installation.  By  the  same  token,  large  values 
for  PEVADl  and  PEVAD2  could  reflect  an  own  ship  sonar  sys¬ 
tem  that  is  not  at  peak  performance  or  not  being  properly 
utilized. 

Similarly,  the  other  parameters  can  be  considered 
to  reflect  the  weapon  and  fire  control  systems  capabili¬ 
ties  of  the  two  ships.  Low  probabilities  of  successful 
snap  shots  are  caused  by  the  fact  that  when  a  snap  shot  is 
fired,  the  fire  control  solution  is  usually  incomplete  if 
there  is  a  solution  at  all.  Also,  the  crew  is  under  some 
stress  due  to  the  fact  that  they  know  a  weapon  is  in  the 
water  heading  toward  them. 

Assumpt ions 

I  he  encounter  envisioned  as  being  applicable  to 
this  model  is  that  of  a  barrier  submarine  versus  a  tran- 
sitor,  although  the  simulation  could  be  modified  so  that 
it  could  be  applied  to  other  types  of  encounters  as  well. 
Some  events  which  might  occur  in  such  an  encounter  have 


been  omitted  because  they  do  not  affect  the  outcome  or 
because  they  are  not  pertinent  to  the  objective  of  the 
simulation.  Encounters  where  the  target  detects  own  ship 
and  attacks  prior  to  own  ship  detecting  the  target  and 
situations  where  own  ship  has  been  trailing  the  target  for 
an  extended  period  of  time  have  not  been  considered. 

Each  submarine  was  assumed  to  have  a  war  load  mix 
of  weapons  on  board  and  it  was  assumed  that  the  maximum 
number  of  weapons  either  ship  could  fire  in  any  one  en¬ 
counter  was  two.  This  weapon  mix  encompasses  all  ranges 
applicable  to  the  encounter.  As  with  the  other  capabili¬ 
ties  considered  for  the  two  ship's  involved,  the  types  of 
weapons  assumed  to  be  carried  can  be  tailored  to  the  type 
of  encounter  desired. 

If  neither  ship  is  sunk,  the  encounter  is  assumed 
to  be  terminated.  It  should  be  realized  that  either  one 
or  both  ships  might  continue  to  prosecute  the  encounter, 
however  it  is  considered  that  the  termination  assumption 
is  not  detrimental  to  the  objective  of  this  study. 

Range  Considerations 

As  was  stated  earlier,  the  values  assigned  to  the 
probabilities  associated  with  each  node  reflect  the  rela¬ 
tive  capabilities  of  the  two  submarines  involved.  It  is 
intuitively  reasonable  that  these  parameters  are  also 
functions  of  the  range  between  the  two  ships  and  the  Lime 
contact  has  been  held.  Accordingly,  as  the  range  de¬ 
creases,  the  parameter  values  will  change.  This  property 
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is  compensated  tin;  Lu  the  model  by  dividing  t  lie  encounter 
into  three  regions.  I  he  selection  and  definition  of  these 
regions  was  somewhat  arbitrary.  Jt  would  be  possible  to 
divide  the  encounter  into  several  more  regions  if  a  more 
detailed  analysis  were  required.  However,  the  three  cho¬ 
sen  are  representative  of  the  profile  of  a  typical  encoun¬ 
ter  and  therefore  describe  the  encounter  sufficiently  for 
the  purpose  of  this  model.  It  should  be  emphasized  that 
the  choice  of  the  regions  is  dependent  on  the  encounter 
under  consideration  so  that  the  physical  areas  contained 
in  each  particular  region  may  vary  between  one  type  of 
encounter  and  another.  In  adapting  this  model  to  a  spe¬ 
cific  study,  selection  and  description  of  regions  must 
be  tailored  to  fit  the  objective  of  that  study. 

In  this  problem,  region  one  describes  that  period 
of  the  encounter  during  which  the  transitor  is  first  de¬ 
tected  and  classified.  A  fire  control  solution  is  ob¬ 
tained  which  would  enable  the  launching  of  a  weapon  with 
a  reasonable  expectation  of  a  successful  kill.  Refine¬ 
ment  of  the  fire  control  solution  coupled  with  decreasing 
range  brings  the  encounter  into  region  two.  As  more  data 
is  collected,  the  solution  is  continually  refined  and  up¬ 
dated  in  this  region.  Region  three  starts  when  the  fire 
control  solution  can  no  longer  be  significantly  improved 
and  it  lasts  down  to  the  minimum  weapon  employment  range. 

I  he  re  are  both  advantages  and  disadvantages  in 
launching  a  weapon  as  early  in  the  encounter  as  possible. 
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An  unalerted  target  is  more  vulnerable  than  an  alerted  one. 
However,  the  later  in  the  encounter  an  attack  is  launched, 
the  more  likely  it  is  that  the  target  counterdetects  the 
attacking  ship  and  presses  an  attack  of  its  own.  At  the 
same  time,  however,  a  premature  attack  that  is  unsuccess¬ 
ful  can  alert  a  target  which  previously  had  no  knowledge 
of  the  presence  of  another  ship.  These  factors  will  be 
investigated  by  analysing  the  results  of  attacks  delivered 
in  each  of  the  three  regions  detailed  above. 

Parameter  Values 

Throughout  the  study,  own  ship  was  assumed  to  de¬ 
tect  the  target  first.  The  parameters  which  were  varied 
in  order  to  examine  the  effect  of  different  capabilities 
are  listed  below. 

PATT  1.  Relatively  constant  over  range 

2.  Decrease  slightly  with  range 

3.  Decrease  proportional  to  range 

PCDET  1.  Capability  comparable  to  own  ship 

2.  Capability  slightly  inferior  to  own  ship 

3.  Completely  inferior  to  own  ship 

PCAT  1.  Well  trained  crew;  Sophisticated  equipment 

2.  Medium  capability 

3.  Poor  capability 

The  values  of  these  parameters,  as  listed  in  Table  I,  are 
based  on  operational  experience  and  are  considered  to  be 
reasonable  representations  of  the  capabilities  considered. 
The  values  for  the  parameters  PSS,  PRAT  and  DPAT  reflect 
the  same  capabilities  as  those  considered  for  PATT.  I  he 
values  for  PCA,  PTAS,  PKVADl  and  PEVAD2  are  derived  f-orr 
the  capabilities  listed  for  PCDET  while  P'l'SS  and  UP  IS 
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LABEL  1 

VALUES  01'  NODAL  PROBABILITIES 


CAPABILITIES 

REGION 

PCDET 

FAIT 

PCA 

PCAT 

PTAS 

PTSS 

1 

I 

.88 

.71 

.30 

.50 

.50 

.18 

1 

2 

.93 

.73 

.50 

.60 

.40 

.23 

1 

3 

.98 

.75 

.70 

.70 

.30 

.28 

2 

I 

.70 

.55 

.30 

.30 

.45 

.13 

2 

2 

.80 

.65 

.50 

.40 

.35 

.18 

2 

3 

.90 

.75 

.70 

.50 

.25 

.23 

3 

1 

.  10 

.25 

.30 

.10 

.30 

.08 

3 

2 

.50 

.50 

.50 

.20 

.20 

.13 

3 

3 

.80 

.75 

.70 

.30 

.10 

.18 

PSS 

PRAT 

PEVADl 

PEVAD2 

DPAT 

UPTS 

1 

1 

.20 

.81 

.10 

.15 

.66 

.20 

1 

2 

.25 

.83 

• 

o 

oo 

.25 

.68 

.25 

1 

3 

.  10 

.85 

.05 

.35 

.70 

.30 

2 

I 

.  15 

.65 

.08 

.10 

.50 

.15 

2 

2 

.20 

. 7  5 

.06 

.20 

.  10 

.20 

2 

'3 

.25 

.85 

.03 

.30 

.70 

.25 

3 

1 

.  10 

.35 

.06 

.05 

.20 

.10 

3 

2 

.  15 

.  60 

.04 

.  15 

.45 

.15 

3 

• 

.20 

.85 

.01 

.25 

.70 

.20 

reflect  those  of  PCAT. 


Measure  of  Effectiveness 

There  are  several  measures  of  effectiveness  which 
could  be  considered  in  the  evaluation  of  the  simulation 
results.  Exchange  ratio  between  the  two  ships  involved  is 
the  measure  of  effectiveness  chosen.  It  is  defined  as  the 
ratio  of  total  target  losses  to  total  own  ship  losses. 
Other  factors,  including  the  four  possible  outcomes  pre¬ 
viously  mentioned,  have  also  been  examined. 
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(HALTER  ill 

RESULTS  01-  THE  SIMULATION 

Qualitatively,  the  results  of  the  computer  simu¬ 
lation  were  reasonable.  These  results,  in  the  form  of 
exchange  ratios  and  outcomes,  are  listed  in  Tables  II 
through  VI.  In  order  to  make  meaningful  use  of  these 
results,  however,  they  must  be  capable  of  quantitative 
analysis.  For  the  determination  of  the  effect  range  and 
elapsed  time  have  on  the  outcome,  contingency  table  tests 
were  made  to  determine  the  dependency  characteristics  be¬ 
tween  regions  and  outcomes.  Then  a  brief  comparison  was 
made  between  the  results  of  the  simulation  and  results 
which  might  be  obtained  by  the  employment  of  the  theory 
of  games  of  timing.  Ln  the  analysis  of  the  effect  of 
different  own  ship  and  target  capabilities  on  the  out¬ 
comes  of  encounters,  a  linear  statistical  model  of  full 
rank  was  used  in  an  attempt  to  determine  exactly  what 
part  is  played  by  each  capability. 

f’onti  ngeney  Tab  I  o  T<  _•  t 

The  tabulat  '<>n  of  exc  har.go  ratios  in  Table  II 
for  given  ship  capabilities  and  regions  leads  to  the  im¬ 
mediate  ((inclusion  that  as  the  elapsed  time  of  the  en¬ 
counter  increases,  the  exchange  ratios  become  Less  favor¬ 
able  to  own  ship.  It  is  intuitively  reasonable  to  assume, 
from  i  h(  -e  tei-ul!-.  that  i!  would  be  advantageous  to  own 
ship  to  preset  ut  <  t  he  at  tail  in  rt  y.  i  on  one  or  as  earl\  in 


TABLE  II 


TABULATION  OF  EXCHANGE  RATIOS 
vs . 

REGIONS  AND  CAPABILITIES 


PCDET 

PCAT 

PATT 

REGION  I 

REGION  II 

REGION 

1 

1 

I 

3.35 

1.32 

.80 

1 

1 

2 

3.30 

1.55 

.87 

1 

1 

3 

2.20 

1.32 

.74 

1 

2 

1 

4.32 

2.25 

1.18 

1 

2 

2 

4.94 

2.46 

1.05 

1 

2 

3 

3.62 

1.93 

1.05 

1 

3 

1 

7.74 

3.54 

1.97 

I 

3 

2 

9.42 

3.52 

1.77 

1 

3 

3 

7.29 

4.18 

1.73 

2 

1 

1 

4.05 

1.94 

1.02 

2 

1 

2 

4.08 

2.29 

.96 

2 

1 

3 

3.35 

1.86 

.92 

2 

2 

1 

5.06 

2.65 

1.29 

2 

2 

2 

5.99 

2.88 

1.45 

2 

2 

3 

6.05 

2.66 

1.28 

2 

3 

1 

7.87 

4.64 

2.10 

2 

3 

2 

10.75 

4.38 

2.09 

2 

3 

3 

8.93 

4.71 

2.51 

3 

I 

1 

9.24 

3.45 

1.19 

3 

1 

2 

10.42 

3.39 

1.22 

3 

1 

3 

9.76 

2.98 

1.25 

3 

2 

1 

9.15 

4.53 

1.45 

3 

2 

2 

11.71 

4.65 

1.54 

3 

2 

3 

11.46 

5.60 

1.59 

3 

3 

1 

8.45 

6.09 

2.46 

3 

3 

2 

12.54 

7.20 

2.72 

3 

3 

3 

16.64 

7.51 

2.77 
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!  MM  I  I.A  1‘  I  ON  uE  slk;<  i-  SSPUl.  ATT  AC. Kb 
vs  . 

RE‘0  1  i.'NS  AN!)  ' ‘,A PA B !  I  \  11ES 


pcdh  r 

PGA  1 

PA  I  1 

RK.IUN  1 

REGION  11 

REGION  ill 

l 

1 

1 

r)b9 

361 

247 

1 

1 

2 

940 

401 

278 

1 

1 

3 

345 

356 

267 

l 

2 

1 

565 

419 

284 

l 

2 

2 

535 

434 

274 

l 

2 

3 

334 

384 

278 

l 

3 

1 

577 

451 

356 

l 

3 

2 

566 

417 

325 

l 

3 

3 

343 

407 

316 

2 

l 

1 

611 

460 

296 

2 

1 

2 

587 

497 

306 

2 

1 

J 

381 

419 

309 

2 

2 

1 

627 

470 

335 

2 

2 

/ 

600 

486 

353 

2 

2 

J 

32  5 

4  35 

349 

2 

J 

1 

653 

483 

362 

2 

3 

2 

612 

510 

369 

2 

3 

3 

390 

464 

39  3 

J 

l 

I 

K  ]  6 

617 

360 

3 

1 

2 

/  b  J 

588 

379 

3 

1 

3 

4  ) 

529 

408 

3 

> 

i 

81  8 

641 

357 

3 

2 

/ 

/4  5 

624 

390 

3 

2 

3 

454 

570 

396 

') 

J 

3 

1 

808 

643 

427 

3 

> 

2 

7  67 

604 

442 

3 

3 

3 

455 

588 

436 

. . . 


TABLE  IV 


TABULATION  OF  BOTH  SHIPS  SUNK 
vs . 

REGIONS  AND  CAPABILITIES 


PCDET 

PCAT 

PATT 

REGION  I 

REGION  II 

REGION  III 

1 

1 

1 

81 

137 

154 

1 

1 

2 

48 

91 

136 

1 

1 

3 

29 

62 

84 

1 

2 

I 

78 

71 

120 

1 

2 

2 

38 

63 

117 

1 

2 

3 

17 

41 

66 

1 

3 

1 

58 

52 

74 

I 

3 

2 

38 

40 

67 

I 

3 

3 

14 

24 

40 

2 

1 

1 

81 

108 

165 

2 

1 

2 

57 

79 

112 

2 

1 

3 

21 

45 

91 

2 

2 

1 

82 

87 

121 

2 

2 

2 

35 

49 

88 

2 

2 

3 

12 

38 

72 

2 

3 

1 

71 

60 

65 

2 

3 

2 

33 

51 

64 

2 

3 

3 

21 

26 

43 

3 

1 

1 

80 

97 

143 

3 

1 

2 

62 

67 

116 

3 

1 

3 

15 

44 

93 

3 

2 

1 

88 

93 

133 

3 

2 

2 

51 

64 

90 

3 

2 

3 

16 

35 

60 

3 

3 

1 

16 

88 

80 

3 

3 

2 

48 

44 

58 

3 

3 

3 

11 

28 

32 

23. 


1  ABLE  V 


TABULA'LL  ON  OF  BOTH  SHIPS  SURVIVE 
vs . 

REGIONS  AND  CAPABILITIES 


PCDET 

PCAT 

pat  r 

REGION  1 

REGION  II 

REGION  III 

1 

1 

I 

237 

2bl 

250 

1 

1 

2 

282 

282 

248 

1 

1 

3 

485 

327 

257 

1 

2 

1 

28b 

363 

374 

1 

2 

2 

349 

364 

354 

1 

2 

3 

569 

396 

395 

1 

3 

1 

341 

407 

426 

1 

3 

2 

381 

453 

454 

1 

3 

3 

608 

490 

478 

2 

1 

I 

218 

247 

251 

2 

1 

2 

255 

251 

257 

2 

1 

3 

499 

331 

257 

2 

2 

1 

233 

320 

311 

2 

2 

2 

294 

328 

342 

2 

2 

3 

561 

387 

321 

2 

3 

1 

255 

400 

435 

2 

3 

2 

328 

362 

424 

2 

3 

3 

564 

432 

433 

3 

1 

L 

87 

176 

219 

3 

I 

2 

160 

219 

215 

3 

I 

3 

466 

279 

192 

3 

2 

1 

83 

197 

305 

3 

2 

2 

187 

228 

299 

3 

2 

3 

505 

322 

318 

3 

3 

I 

85 

237 

367 

3 

3 

2 

168 

306 

374 

3 

3 

3 

517 

330 

395 

24. 
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TABLE  VI 


TABULATION  OF  OWN  SHIP  SUNK,  TARGET  SURVIVES 

vs . 

REGIONS  AND  CAPABILITIES 


PCDET 

PCAT 

PATT 

REGION  I 

REGION  II 

REGION  III 

1 

I 

1 

113 

241 

349 

1 

1 

2 

130 

226 

338 

1 

1 

3 

141 

255 

392 

1 

2 

1 

71 

147 

222 

1 

2 

2 

78 

139 

255 

1 

2 

3 

80 

179 

261 

1 

3 

1 

24 

90 

144 

1 

3 

2 

25 

90 

154 

1 

3 

3 

35 

79 

166 

2 

1 

1 

90 

185 

288 

2 

1 

2 

101 

173 

325 

2 

1 

3 

99 

205 

343 

2 

2 

1 

58 

123 

233 

2 

2 

2 

71 

137 

217 

2 

2 

3 

52 

140 

258 

2 

3 

1 

21 

57 

138 

2 

3 

2 

27 

77 

143 

2 

3 

3 

25 

78 

131 

3 

1 

1 

17 

110 

278 

3 

1 

2 

17 

126 

290 

3 

1 

3 

36 

148 

307 

3 

2 

1 

11 

69 

205 

3 

2 

2 

17 

84 

221 

3 

2 

3 

25 

73 

226 

3 

3 

1 

11 

32 

126 

3 

3 

2 

17 

4b 

126 

3 

3 

3 

17 

54 

137 
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the  encounter  a<  j c s s  •  t)  i  e  .  In  order  to  support  this  con¬ 
clusion,  it  would  be  desirable  to  show  that  the  occurrence 
of  the  possible  outcomes  of  the  encounters  are  dependent 
on  the  region  in  which  the  attack  is  launched.  Tests  of 
independence  in  contingency  tables  as  discussed  in  refer¬ 
ence  4  provide  a  convenient  means  of  examining  these  de¬ 
pendency  properties.  1 f  an  encounter  was  to  be  considered 
an  experiment  whose  results  could  be  classified  by  two 
attributes,  namely  region  and  outcome,  then  a  contingency 
table  would  be  an  appropriate  model.  Such  a  model  is  il¬ 
lustrated  below. 

SUCCESSFUL  BOTH  SHIPS  BOTH  SHIPS  OWN  SHIP 


ATTACKS 

SUNK 

SURVIVE 

SUNK 

REGION 

I 

X11 

X12 

X14 

X14 

REGION 

II 

X2J 

X22 

X2  3 

X24 

REGION 

III 

X31 

X32 

X33 

X34 

The  X^'s  represent  the  number  of  appropriate  outcomes  of 
the  encounter  while  the  random  variable 

4  3  jx..  -  n(X../n)(X. ./n)J  2 


z  r 

i=i  i=i 


n(x../n)(X. ./n) 


X^.  =  Sum  of  the  elements  of  row  i. 

X.j  =  Sum  of  the  elements  of  column  j. 
n  =  Number  of  observations  (encounters), 
has  an  approximate  (hi -square  distribution  with  five  de¬ 
grees  of  freedom.  If  the  value  of  the  above  expression  is 
greater  than  the  value  in  the  Chi-square  table  for  the  de¬ 
sired  significance  level  then  the  hypothesis  that  the  out- 
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comes  are  independent  of  the  regions  is  rejected  at  that 
significance  level.  Using  the  results  of  three  thousand 
encounters  for  each  of  the  twenty-seven  different  combi¬ 
nations  of  the  capabilities,  it  was  found  that  the  hypo¬ 
thesis  of  independence  was  rejected  at  any  level  of  signif¬ 
icance  tabulated.  The  conclusion  is  therefore  drawn  that 
the  independence  assumption  should  not  be  accepted  and 
that  the  intuitive  analysis  of  the  table  of  exchange  ratios 
is  valid. 

Games  of  Timing 

The  proposition  that  it  is  desirable  to  shoot  as 
early  in  an  encounter  as  possible  is  also  borne  out  by  the 
theory  of  games  of  timing,  or  duels,  as  discussed  in  ref¬ 
erence  1.  The  type  of  encounter  which  is  the  subject  of 
this  study  would  be  classified  as  a  noisy  duel  with  many 
bullets.  In  such  a  duel,  the  antagonists  should  shoot  as 
soon  as  the  following  equation  is  satisfied: 

P(t)  - - - - 

m( t)+n( t) 

P(t)  is  the  accuracy  function  or  probability  of  successful 
attack,  while  ro(t)  and  n(t)  are  the  number  of  weapons 
available  on  each  ship  for  the  particular  encounter.  Al¬ 
though  this  particular  model  assumes  equal  accuracies  for 
both  ships,  it  does  provide  a  means  for  studying  an  encoun¬ 
ter  of  this  type.  Assuming  that  each  ship  has  at  least  two 
weapons,  the  value  of  P(t)  which  satisfies  the  equation  is 
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less  than  or  equal  to  .25.  This  value  occurs  in  region 
one  for  all  of  the  encounters  which  were  considered.  An 
interesting  study  would  be  the  application  of  the  theory 
of  stochastic  duels,  extended  to  include  varying  single 
shot  kill  probabilities  as  developed  in  reference  2,  to 
this  submarine  versus  submarine  situation. 

Linear  Statistical  Model 

In  order  to  study  the  effect  that  different  own 
ship  and  target  capabilities  have  on  the  outcome  of  an 
encounter,  use  was  made  of  a  linear  statistical  mod^l  of 
full  rank.  In  this  type  of  model,  there  is  a  known  out¬ 
come  assumed  to  be  a  random  variable  which  is  a  function 
of  p  unknown  parameters.  If  there  are  n  outcomes  then  the 
known  quantities  form  a  matrix  and  the  unknown  parameters 
form  a  vector.  Written  in  matrix  notation,  the  linear 
model  becomes: 

Y  =  XB  +  e 

Y  is  an  nxl  vector  of  known  outcomes. 

X  is  an  nxp  matrix  of  known  quantities.  The  rank 
X  is  p . 

B  is  a  pxl  vector  of  unknown  parameters. 

e  is  a  pxl  vector  of  error  terms. 

A  least  squares  solution  to  this  equation  yields  a  best 
linear  unbiased  estimator  of  the  vector  B.  Written  in 
matrix  form,  this  estimator,  B,  is: 

B  =  (X’X)"1X'Y 
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Returning  to  the  objective  of  this  study,  it  is  obvious 
that  if  the  outcomes  of  the  encounters  and  the  capabili¬ 
ties  of  the  ships  involved  could  be  fitted  into  such  a 
model,  then  the  values  for  B  should  be  good  estimates  of 
the  relative  effects  which  the  various  capabilities  have 
on  the  outcomes. 

Initial  attempts  at  this  method  were  unsuccessful. 
Five  Y  vectors,  dimensioned  16x1,  were  examined  using  X 
matrices  dimensioned  16x8  and  16x9.  The  five  Y  vectors 
corresponded  to  the  possible  outcomes  of  an  encounter  and 
the  exchange  ratios  as  detailed  below: 

Y^  -  Successful  Attacks 

Y2  -  Both  Ships  Sunk 

Y^  -  Both  Ships  Survive 

Y^  -  Own  Ship  Sunk,  Target  Survives 

Y^  -  Exchange  Ratios 

The  values  of  the  components  of  the  vectors  Y^  through 
Y^  were  divided  by  1000  to  give  a  relative  frequency  of 
occurrence  for  any  given  set  of  capabilities.  Eighteen 
encounters  were  selected  out  of  the  eighty-one  possible 
combinations  of  capabilities  and  regions.  The  encounters 
chosen  included  six  different  capability  mixes  with  at¬ 
tacks  occurring  in  each  of  the  three  regions.  The  di¬ 
mension  was  reduced  from  eighteen  to  sixteen  to  make  the 
X  matrix  nonsingular.  The  eight  column  vectors  of  the  X 
matrix  corresponded  to  the  appropriate  values  for  PCDET, 
PATT,  PC AT ,  PSS,  PEVAD1  and  PEVAD2 .  The  other  nodal  values 
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were  not  included  because  they  formed  linear  combinations 
of  the  other  columns  and  would  have  resulted  in  a  singular 
matrix.  In  order  to  reduce  the  restrictions  on  the  solu¬ 
tions  a  unit  column  vector  was  added  to  the  X  matrix  bring- 

A 

ing  its  total  dimension  to  16x9.  The  values  for  B  arrived 
at  by  this  computation  were  not  satisfactory.  They  did 
not  reflect  the  relationship  which  should  exist  between 
the  parameters  and  the  outcomes.  This  was  borne  out  by 
the  fact  that  when  the  matrix  multiplication,  XB,  was  per¬ 
formed,  it  did  not  yield  an  estimate  of  the  Y  vector  which 
could  be  correlated  in  any  manner  with  the  actual  compo¬ 
nents  of  Y.  The  many  unsuccessful  attempts  to  fit  the 
parameter  and  outcome  values  to  the  model  indicated  a  con¬ 
clusion  that  the  factors  affecting  the  outcome  of  an  en¬ 
counter  do  not  have  a  linear  relationship  to  that  outcome. 

However,  there  are  three  factors  which  were  pre¬ 
viously  assumed  to  embody  the  capabilities  reflected  by  the 
other  parameter  values.  These  three  factors  are  PCDET, 

PATT  and  PCAT.  Column  vectors  of  these  values  plus  a  unit 
column  vector  were  used  to  form  a  new  X  matrix.  The  ele¬ 
ment  of  B  corresponding  to  the  unit  column  vector  could  be 
considered  to  incorporate  the  terms  deleted  from  the  origi¬ 
nal  X  matrix  so  that  the  value  reflects  the  cumulative 
effects  of  the  deleted  terms.  Since  the  previous  investi¬ 
gation  into  the  effect  that  range  has  on  the  outcome  of  an 
encounter  indicated  that  the  optimum  time  to  launch  an  at¬ 
tack  was  in  region  one,  only  that  region  was  considered  in 


30. 


the  new  analysis.  The  linear  model  was  then  written: 

Y  =  Bj^l  +  B2X2  +  B3X3  +  B4X4  +  e 
B2 ,  B3 ,  and  B4  are  the  elements  of  the  B  vector.  The 
unit  column  vector  is  represented  by  1,  X2  is  the  column 
vector  of  PCDET  values,  X3  is  the  column  vector  of  PATT 
values  and  X4  is  the  column  vector  of  PCAT  values.  The 
values  for  B  resulting  from  these  computations  are  shown 
in  Table  VI I . 


TABLE  VII 

VALUES  OF  B  FOR  THE  DIFFERENT  Y  VECTORS 


*1 

Y1 

.50281 

Y2 

.01544 

Y3 

. 44448 

Y4 

.03727 

Y5 

8.36105 

^2 

-.30239 

-.02100 

.24811 

.07527 

-4.96775 

^3 

.47109 

.10424 

-.47537 

-.09996 

3.38144 

.00874 

.00562 

-.14913 

.13477 

.5.22542 

The  estimates  of  the  Y  vectors  resulting  from  the  multipli- 

A 

cation  XB  coincided  more  closely  with  the  actual  values  of 
Y  than  did  the  values  obtained  with  the  larger  X  matrix.  In 
addition,  nine  sets  of  capabilities  which  were  not  used  in 
the  determination  of  ^  were  substituted  into  the  linear  mo¬ 
del.  The  estimates  of  Y  coincided  extremely  well  with  the 
actual  values  of  Y  obtained  from  the  simulation.  The  results 
of  these  computations  are  contained  in  Tables  VIII  and  IX. 
This  should  indicate  that  a  linear  relationship  does  exist 
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TABLK  Vi  1 1 


RESULTS  OF  THE  LINEAR  STATISTICAL  MODEL 


Successful  Attacks 


Y' 

=  ( 

.569 

.540 

.345 

.611 

.816 

.565 

.577 

.627 

.808) 

Y' 

=  ( 

.576 

.500 

.359 

.630 

.811 

.574 

.572 

.628 

.808) 

1y' 

-Y' 

1  -  ( 

.007 

.040 

.014 

.019 

.005 

.009 

.005 

.001 

.000) 

Both 

Shi^s  Sunk 

Y' 

=  ( 

.081 

.048 

.029 

.081 

.080 

.078 

.058 

.082 

.096) 

-  ( 

.074 

.057 

.026 

.078 

.090 

.073 

.072 

.076 

.088) 

|y' 

-Y ' 

1  =  < 

.007 

.009 

.003 

.003 

.010 

.005 

.014 

.006 

.008) 

Both 

Ships  Survive 

Y' 

=  ( 

;.237 

.282 

.485 

.218 

.087 

.286 

.341 

.233 

.085) 

Y' 

*  ( 

;  .251 

.327 

.469 

.206 

.057 

.281 

.310 

.236 

.117) 

(y' 

-Y ' 

1=  < 

:.oi4 

.045 

.016 

.012 

.030 

.005 

.031 

.003 

.032) 

Own  Ship  Sunk, 

Target  Survive 

s 

Y' 

=  ( 

.113 

.130 

.141 

.090 

.107 

.071 

.024 

.058 

.011) 

Y' 

=  ( 

.100 

.116 

.146 

.086 

.041 

.073 

.046 

.059 

.013) 

|y' 

A 

-Y’ 

1  -  ( 

.013 

.014 

.005 

.004 

.024 

.002 

.022 

.001 

.024) 

Exchange 

Ratio 

s 

Y' 

-  ( 

3.35 

3.30 

2.20 

4.05 

9.24 

4.32 

7.74 

5.06 

8.45 

Y’ 

-  ( 

3.78 

3.24 

2.22 

4.67 

7.65 

4.82 

5.87 

5.72 

9.74 

|v 

-Y' 

1=  < 

.43 

.06 

.02 

.62 

1.59 

.50 

1.87 

.66 

1.29 
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TABLE  IX 


COMPARISON  OF  OTHER  VALUES  OF  Y 


Successful  Attacks 


Y' 

3 

(.535 

.334 

.566 

.587 

.381 

.600 

.375 

.653 

.818) 

as 

(.498 

.357 

.497 

.555 

.413 

.553 

.412 

.626 

.810) 

|Y'-Y'| 

3 

(.037 

.023 

.069 

.032 

.032 

.047 

.037 

.027 

.008) 

Both  ; 

Ships  Sunk 

Y' 

* 

(.038 

.017 

.038 

.057 

.021 

.035 

.012 

.071 

.088) 

Y' 

s 

(.056 

.025 

.055 

.061 

.030 

.060 

.028 

.075 

.089) 

lY'-Y'l 

s 

(.018 

.008 

.017 

.004 

.009 

.025 

.016 

.004 

.001) 

Both  Ships  Survive 

Y' 

3 

(.349 

.569 

.381 

.255 

.499 

.294 

.561 

.255 

.083) 

(.357 

.499 

.386 

.282 

.424 

.312 

.455 

.266 

.087) 

lY'-Y'l 

3 

(.008 

.070 

.005 

.027 

.074 

.018 

.106 

.011 

.004) 

Own  Ship 

Sunk , 

Target  Survives 

Y' 

* 

(.078 

.080 

.025 

.101 

.099 

.071 

.052 

.021 

.011) 

</• 

= 

(.089 

.119 

.062 

.102 

.132 

.075 

.105 

.032 

.014) 

|Y'-Y'| 

a 

(.011 

.039 

.037 

.001 

.033 

.004 

.053 

.011 

.003) 

Exchange 

Ratios 

Y' 

3 

(4.94 

3.62 

9.42 

4.08 

3.35 

5.99 

6.05 

7.87 

9.15) 

A 

Y' 

= 

(4.28 

3.27 

5.33 

4.13 

3.12 

5.18 

4.16 

6.76 

8.70) 

|Y'-Y'| 

3 

(  .  66 

.35 

4.09 

.05 

.23 

.81 

1.89 

1.10 

.45) 
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between  outcomes  and  the  factors  affecting  them,  and,  if 
that  relationship  could  be  found,  a  meaningful  analysis 
of  the  factors  involved  could  be  made. 

The  calculations  associated  with  the  linear  sta¬ 
tistical  model  were  made  on  an  IBM  360  Series  computer. 
The  program  for  those  calculations  is  contained  in  Appen¬ 
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CHAPTER  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions 

This  study  has  attempted  to  model  an  encounter 
between  two  submarines  and  to  quantitatively  analyse  the 
factors  involved.  The  logic  tree  and  the  computer  simu¬ 
lation  together  form  a  representation  of  encounters  which 
can  be  made  applicable  to  a  variety  of  own  ship  and  tar¬ 
get  capabilities.  A  large  sample  of  encounters  can  thus 
be  analysed. 

The  results  of  the  model  indicate  the  desirability 
of  prosecuting  an  attack  in  region  one  instead  of  regions 
two  or  three.  The  demonstration  of  the  dependency  of  the 
outcomes  on  the  region  in  which  the  attack  is  launched 
lend  added  significance  to  this  conclusion,  thus  support¬ 
ing  present  tactical  doctrine.  This  support  also  tends 
to  add  credence  to  the  validity  of  the  model. 

The  use  of  the  linear  statistical  model  and  the 
results  obtained  from  it  were  disappointing  at  first.  How¬ 
ever,  a  linear  unbiased  estimator  which  could  accurately 
estimate  the  original  Y  vector  was  obtained  using  a  re¬ 
duced  size  X  matrix.  This  final  estimator  did  reflect 
reasonable  representations  of  the  effects  the  different 
factors  should  have  on  the  outcomes  of  the  encounters. 
Further  refinement  of  this  technique  should  yield  a  method 
for  obtaining  a  meaningful  analysis  of  the  effect  proposed 
systems  would  have  on  ship  capabilities. 
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Recommendations 


The  model  proposed  in  this  thesis  is  not  intended 
to  replace  other  types  of  models  now  in  use,  nor  is  it 
designed  to  solve  all  problems.  However,  for  problems  in¬ 
volving  two  opposing  forces  and  whose  characteristics  can 
be  applicably  quantified,  this  basic  model  presents  a  meth¬ 
odology  which  could  be  an  effective  tool  in  the  hands  of 
an  analyst. 

Application  of  this  model  to  an  actual  situation 
immediately  poses  the  question  as  to  where  values  for  the 
nodal  probabilities  can  be  obtained  or  in  what  manner  can 
they  be  determined?  The  simulation  and  the  analysis  of 
its  results  become  an  academic  exercise  if  they  cannot  be 
compared  with  actual  results  and  applied  to  actual  prob¬ 
lems.  Fortunately,  there  is  a  wealth  of  information  con¬ 
tained  in  the  reports  of  exercises  conducted  by  various 
submarine  commands  and  by  SUBMARINE  DEVELOPMENT  GROUP  TWO. 
Since  the  subject  of  this  thesis  was  to  propose  and  devel¬ 
op  a  model,  actual  exercise  data  was  not  used,  although 
the  specific  values  desired  could  be  extracted  from  such 
exercise  data  and  therby  apply  the  model  to  a  "real  world" 
situation. 

This  study  has  pointed  to  several  areas  where  fur¬ 
ther  development  on  the  thesis  level  should  prove  to  be 
fruitful.  These  areas  are  detailed  below: 

1.  This  model  should  be  applicable  to  a  variety  of 
types  of  encounters.  One  that  immediately  comes 
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to  mind  is  an  aircraft  versus  aircraft  encounter. 

A  further  extension  would  be  an  attempt  to  apply 
it  to  encounters  between  opposing  forces  where  each 
force  is  composed  of  multiple  units. 

2.  Further  development  of  the  linear  statistical  model 
is  needed  to  obtain  a  best  linear  unbiased  estima¬ 
tor  which  more  closely  reflects  the  effects  all  of 
the  factors  have  in  determining  the  outcome  of  an 
encounter . 

3.  According  to  the  Central  Limit  Theorem,  the  large 
number  of  observations  enable  the  assumption  to  be 
made  that  the  e^'s  of  the  linear  statistical  model 
are  normally  distributed  random  variables.  Under 

this  assumption,  hypothesis  tests  on  and  confidence 

/\ 

intervals  for  the  values  of  B  can  be  obtained. 

Such  statistical  tests  should  be  made  in  order  to 
gain  more  insight  into  the  problem. 

4.  An  attempt  to  apply  the  theory  of  stochastic  duels 
with  varying  single  shot  kill  probability,  as  de¬ 
tailed  in  reference  2,  to  the  submarine  versus  sub¬ 
marine  encounter  should  also  prove  to  be  profitable. 
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CASE  II  |  CASE  1 


n 


s 


ID?sIP8+ 1 


*P£MD 


NMuPT: 


IW7cIW7+\ 

JdI3*IDI3+| 


APPENDIX  C 


COMPUTER  PROGRAM  FOR  l HE  SIMULATION  OF  A 
SUBMARINE  VERSUS  SUBMARINE  ENCOUNTER 


DIMENSION  PCOEI  I  3, 3 ) , PATT ( 3 , 31 , PCA ( 3, 3 1 » PC AT ( 3 , 3 ) , 

1 P  T  A  S  ( 3 . 3 )  ,PTSS(3.3).PSS(3.3»,PPAT(  3.3)  , PEVAD1 ( 3, 3) , 
1PEVAP2,3 ,3}  ,OPAT t  3,3 ) .UPT$( 3,3) ,  IO< 15 ) , IW<  15 ) , 
IIL(15),PNI?5),ERI9,3,3) 

0UM=U*M<  01 

I K  =0 
JK=0 


R  F  A  0  (  c  !  i 

(  ir-rijETI  I,  J) 

,  J=1 

,3), 

1=1,  3) 

RE  AO  I 5,Z  ) 

( ( PATT 

(  I  ,  J) 

,J=1 

,3), 

1  =  1,3) 

REAL)  (  5 . 1  ) 

I  (  OCA 

(  I,  J) 

f  J  =  1 

♦  3  )  , 

1=1,3) 

Rr  AO (5,4  » 

I  (  PC  A  T 

I  I,  J) 

,  J=1 

,3)  , 

1  =  1,3) 

RE  AD ( 5 , 5  ) 

I ( FT  AS 

(  I  ,  J) 

f  J—  1 

.3), 

1=1,3) 

REAQI^C  ) 

!  •  PTSS 

(  I  ,  J) 

,J  =  1 

i  3  )  , 

1  =  1,3) 

BFAr'5,7  ) 

I  *  PSS 

’  I  ,  J  ) 

,  )  =  1 

,3), 

1=1 , 3» 

1  =  1.3) 

READ<5,6  I 

( ( HKA1 

(I  ,  J) 

,  J=  1 

,3), 

REA'M%S  ) 

«  !  OFVAOl II  ,  J 

)  ,  J  = 

1,3) 

•  I  = ! ,3) 

ocun f  F . i n i 

(  I PfcvAD2 (  I  ,  J 

I  *  »  = 

1,3) 

,I’i, 3) 

RE  ADI  5, i 1 1 

I  I jFm 1 

i  I  ,  J 

)  ,  J  = 

1,3) 

,1=1.3) 

OF  AO ( 5 , 1 2 ) 

p;;rm,m(cft 

( (OPTS 

(  I  ,  J 

)  ,  J  = 

1,3) 

,1=1,3) 

1 

.  ?) 

2 

F0Rm,-,t<9F 

.2» 

3 

FORMAT ' 9 1  ' 

.2! 

4 

F0RMAM9F3 

.2) 

5 

FGRMATtQP3 

.  >) 

6 

FORMATf QF3 

.2) 

7 

FORMAT (9F3 

.2) 

^CRmaTI^F^ 

.?: 

9 

F0R**ATi9C: 

.2* 

1  0 

FO  RM  AT  1 QF  j 

.2) 

11 

FORMAT ( RF3 

.2) 

12 

FORMAT  I SF3 
WR I TE ( 6 , 13 

.1) 

\ 

13 

FORMAT! 1H1 ) 

L  IS  THF  INDEX  for  VARIATION  OF  PCDET  CAPABILITY. 

K  IS  THE  INDEX  FOR  VARIATION  OF  9C,  AT  CAPABILITY. 

I  IS  THE  !“tiFX  FOR  VARIATION  OF  PATT  CAPABILITY. 

J  IS  ThC  H'OEX  FCr  VARIATION  OF  R 1 3  T  ON  IN  WHICH 

ATTACK  IS  LAUNCHED. 

DO  5000  L  =  i » 3 
DO  40CC  K-) ,3 
JK=JK+I 
DO  3000  1*1,3 
DO  2CC0  J= 1  *  3 
DC  20  11*1,15 
10  (  I  1 )  =0 
I L  (  I  1  )  -G 

20  IW(ll)=C 
IOT  =  0 

I L  T  =  0 
Ik7=o 
l  DT 1 *0 
1 0  T  2  =  0 
IL  T1 =0 
I L  T  2  *  0 
I WT 1 =G 
1WT2-0 

:cn  =  o 

on  90  ji=i,iooo 

GENERATION  OF  25  RANDOM  NUMBERS  FOR  FACH  ENCOUNTFP 
DO  21  11=1,25 
RNI I  1 )=URNi 1 ) 

21  CONTINUE 
ENCOUNTER  STARTS 

IF(RNIl) .LE.PCDETIL, J)  > GO  TO  5C 
CASF  I 


-+  )  . 


IF(RN(2).LF.PATT(I,J) 1G0  TO  22 
IF(RN<5)  .LE.PTASU  1GO  TO  24 
IF(RN(6) .LE.PEVAOll I .  J ) 1  GO  TO  27 
IF(RN(7).LE.PRAT(I,jhGO  TO  28 
10(51=10(51+1 
GO  TO  <90 

22  IF(RN( 3) .LE.PTASl I ,J) 1GO  TO  23 
iwm=iwm  +  i 

GO  TO  90 

23  IF(RN(4) .LE.PTSS( I » J ) )G0  TO  31 
IW<  2 1=1 W(2 1+1 

GO  TO  90 

24  IF f RN( 8 1  .LE.PTSSU  1G0  TO  25 
IF(RN( 10).LE.PRAT( I, J  ))  GO  TO  26 
10(31=10(31+1 

GO  TO  90 

25  IF(RN(9) .  LE . PRAT ( I ,  J 1 )G0  to  32 
IL(ll=IL(i)+l 

GO  TO  90 

26  IW(3)=IW(3)+1 
GO  TO  90 

27  10(41=10(41+1 
GO  TO  90 

28  IW(4)=IW(4)+l 
GO  TO  90 

31  10(11*10(11+1 
GO  TO  90 

32  10(21=10(21+1 
GO  TO  90 

C  CASE  II 

50  I CD= ICO+ 1 

I F ( R N ( 11 ) ,LE.PCA( I  •  J  1  )G0  TO  70 
IF(RN(24).LE.PEVA02( I . J) )G0  TO  62 
IF(RN(15).IE»0PAT( I » J  1 1 60  TO  58 
IF(RN(18).IE.PTAS< I . J 1 1  GO  TO  54 
IF(RN(19). LE.PEVAOll I •  J  1 ) GO  TO  52 
IF(RN(20).LE.PRAT( I, J))GO  TO  53 
I0( 101=1 D( 101+1 
GO  TO  90 

52  I0(9|*I0(9)+1 
GO  TO  90 

53  IW(91=IW(9)+1 
GO  TO  90 

54  IF ( RN( 21 J.IE.UPTS( I, J I) GO  TO  55 
IF(RN(221.LE.PRAT(  I,  JMGO  TO  57 

io(in*ro(ii)  +  i 

GO  TO  90 

55  IF(RN(23).LE.PRAT(  I»J) ) GO  TO  56 
IL(3)=IL(3)+l 

GO  TO  90 

56  I0(12)=I0( 121+1 
GO  TO  90 

57  IW(8)=IW(8)+1 
GO  TO  90 

58  IF(RN( t61.LE.PTAS( I, J1IG0  TO  5« 
IW(6)=IW(6)+1 

GO  TO  90 

59  IF ( RN( 1 7 1 . LE .UPTS(  I » J 1 1  GO  TO  60 
lwm«lM(7)  +  l 

GO  TO  90 

60  10(131*10(131+1 
GO  TO  90 

62  10(81*10(81+1 
GO  TO  90 

70  IF  ( RN(  1 2  I  •(.  E  •  PC  AT  (  K  t  J  1  )  GO  TO  71 
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IF  IRISH  13)  .LE  .PSS  (  !  ,J)  )GD  TO  73 
104  7  >  =  ID*  7 ) ♦  1 

go  in  <»o 

71  IF ( RN 4 14J.LE.PSSl I , J) ) GO  TO  72 
IL(2)  =  IL(;:)  +  1 

GO  TO  90 

72  ID(6)*ID(6)+1 
GO  TO  90 

73  IW(5)*IW(5)+1 

90  CONTINUE 

00  91  K1 =1 ,5 

91  IDT1*I0T1+I0(K1  ) 

DO  92  Kl*6 , 1 3 

92  !DT2*IQT2*iO(K1 ) 

DO  93  K 1 =1  •  4 

93  IWTl=IWTi+f WIKI) 

00  94  K1 =5 . 9 

94  IMT2=IMT2«-IM(K1) 

ILT1=IL( 1) 

ILT2*IL(2)+IL( 3) 

IDL1  *  I  D(  1 )  ♦  I  D(  2  ) 

I0L2*  ID( 6 1 ♦ I  04 12)  +  ID< 131 

IDM2  = I D( 7)+ID(8)+ID{9)+ID( 1 0 ) ♦ I  0  4 11) 

IDW1-IDI 3)MU(4)  +  ID(5) 

IDT* I  DTI ♦ I0T2 

ILT*ILTUILT2 
I WT= I WT1 ♦ I WT  ? 

I  SUM* ID* 1)  +  ID(2)«-ID(6)  +  ID( 1 2 )  + 1 D  ( 13) 
SUMl= I WT+I  S-J* 

SUM2*lLTfISU^ 

EXR*  ( SUM 1 1 /  (  SUM2 ) 

IDL*ISUM 
IDM* IDT- I  SUM 
ER ( JK • I , J) *E  XR 
MRITE(6,200) J. I iK,L 
MR ITE ( 6*  2C1 ) 

WRITE  <  6*202  I  I H( 1 ) 

MR  I TE ( 6,203 ) IW(2) 

MRITE (6,204) I W( 3 ) 

MR  I TE I  6,205 ) I W( 4 ) 

MR  I TE I  6 »  206 ) I D (  1  ) 

MR  I TE ( 6  * 207 )  ID*  2  I 
MR!TE*6,208) I D ( 3 i 
MR  I TE ( 6 ,  209 )  I D  *4  I 
MR  I TE ( 6*  210 ) I D *  5  * 

MRITE (6,211 ) II ( 1  ) 

MRITE (6*212 ) 

MR  I TE 1 6 , 213 1  I W ( 5) 

MRITE<6,214) I M ( 6 ) 

MR  I TE I  6 , 215 ) I M ( 7  I 
MR  I TE  *  6, 216 ) IW<8> 

MRITE  6,217 ) IM(9 ) 

HRITE*6, 218)10(6) 

WRITE ( 6,219) ID( 7 ) 

MR  I TE  *  6 , 220  > 10(8) 

MRITEI6, 221) I D ( 9 ) 

MRITE(6, 222)10(10) 

MR  I te  4  6 , 223  ) I D(  1 1 ) 

MRITE(6,224) ID( 12) 

MR  I TE  4  6 , 225  ) I D ( 13) 

MRITE  6,226) 11(21 
MR  I TE ( 6, 227 ) IL ( 3  ) 

MR  I TE ( 6 , 2 34 ) I WT 
MR  I TE  4  6 , 235  ) IDL 
MR  I TE ( 6 , 235 1  ) I OM 


taR I TE ( 6 , 236) ILT 
WRITE(6,237) ICO 
WRITE(6*23R)  EXR 
WRITE(6,239) 

20C  FORMAT(15X,6HREGION, I  3 , 5X , 4HP ATT  ,  I  3 , 5X , 4HPC AT , I  3,  5X, 
15HPCOET, 13/) 

201  FORMAT(15X*6HCA$E  I./) 

202  FORMAT(15X»3 OH SUCCESSFUL  ATTACKS  (I)  =,I4/) 

203  FORMAT(15X*30H SUCCESSFUL  ATTACKS  ( 2 )  =,I4/> 

204  F0RMAT(15X, 30H SUCCESSFUL  ATTACKS  (3)  =,I4/) 

205  FORWAT(15X,30HSUCCESSFUL  ATTACKS  (4)  =,I4/) 

206  FORMAT(15X,30HBOTH  SHIPS  SUNK  (I)  =  s 1 4/ ) 

207  FORHAT(15X,30HBOTH  SHIPS  SUNK  (2)  * ,  14/) 

208  FORHAT(15X,30HBOTH  SHIPS  SURVIVE  (1)  = , 14/) 

209  F0RMAT(15X,30HB0TH  SHIPS  SURVIVE  (4)  = ♦ 1 4/ ) 

210  FORMAT(15X,30HBOTH  SHIPS  SURVIVE  (5)  =,l4/> 

211  FORMAT! 15X, 30H0WN  SHIP  SUNK  (1)  =  *  14/) 

212  FORMAT(15X,7HCASF  II./) 

213  FORMAT! 15X»  30H SUCCESSFUL  ATTACKS  (5)  =,I4/) 

214  F0RMAT(15X,30HSUCCESSFUL  ATTACKS  (6)  =,I4/) 

215  F0RMAT(15X, 3 OH  SUCCESSFUL  ATTACKS  (7)  =,I4/) 

216  FORMAT! 15X» 3 OH SUCCESSFUL  ATTACKS  (8)  =  *  I  4/ > 

217  FORMAT!15Xf3 OH  SUCCESS  FUL  ATTACKS  (9)  = ,  14/) 

218  FORMAT!15X,30HBOTH  SHIPS  SUNK  (6)  =,I4/) 

219  FORMAT! 15X, 30HB0TH  SHIPS  SURVIVE  ( 7 )  =,I4/) 

22C  FORMAT(15X,30HBOTH  SHIPS  SURVIVE  «8)  =,I4/) 

221  F0RMAT(15X, 30HBOTH  SHIPS  SURVIVE  (9)  =,I4/) 

222  FORMAT(15X»30HBOTH  SHIPS  SURVIVE  (10)  =,I4/) 

223  F0RMAT(15X.30H8OTH  SHIPS  SURVIVE  (11)  *,I4/) 

224  FORMAT(15X,30HBOTH  SHIPS  SUNK  (12)  =,I4/) 

225  FORMAT! 15X.30HBOTH  SHIPS  SUNK  (13)  =,I4/) 

226  F0RMAT(15X, 30HOWN  SHIP  SUNK  (2)  =,I4/) 

227  F0RMAT(15X»  3 OH OWN  SHIP  SUNK  (3)  * . 14/) 

234  F0RMAT(15X»4 OH TOTAL  SUCCESSFUL  ATTACKS 

l  = ,14 / ) 

235  F0RMAT(15X,40HT6taL  BOTH  SHIPS  SUNK 

1  = » 14/ ) 

2351  F0RMATI15X, 40HTOTAL  BOTH  SHIPS  SURVIVE 
1  =,I4/) 

236  FORMAT(15X,40HTOTAL  OWN  SHIP  LOSS, TARGET  SURVIVES 

1  = ,  I 4/  ) 

237  FORMAT! 15X.40HTOTAL  OWN  SHIP  DETECTED  BEFORE  ATTACK 

1  =  , 1 4/ ) 

238  FORMAT(15X.17H EXCHANGE  RATIO  =  , F3 . 2 , 2X , 4HT0  1) 

239  FORMAT ( 1H1 ) 

200C  CONTINUE 
3000  CONTINUE 
400C  CONTINUE 
500C  CONTINUE 

WRITE! 6,600 1 ) 

6001  FORMAT! 1H1 , 43X , 1 5HEXCHANGE  RATIOS//) 

DO  6008  L= 1 , 3 

DO  6006  K=  1 , 3 
IK*I K+l 

WRITE(6,6003)L,K 
WRITE! 6,6004 ) 

DO  6002  1=1,3 

WRITE (6,6005 )I,(ER(IK,I,J),J=i,3) 

6002  CONTINUE 

6003  FORMAT(40X,5HPCDET ,13, 1CX,4HPCAT • I  3/) 

6004  FORMAT! 1 5X, 6HREG ION. 13X  ,  1H1 , 13X , 1H2, 1 3X, 1H3, / ) 

6005  F0RMAT(15X,4HPATT,I3,9X,F5.2,9X,F5.2,9X,F5.2,10X, 


14HT0  1/) 
6006  CONTINUE 
6008  CONTINUE 


48. 


APPEND  LX  D 


sample:  output  of  simulation 

REGION  1  PAI T  1  PCAT  1 


CASE  1 


SUCCLSSLUl.  ATI  APRS 

(1) 

=  42 

SUCCESSFUL  ATTACKS 

(2) 

ii 

cr 

SUCCESSFUL  ATTACKS 

(3) 

=  11 

SUCCESSFUL  ATTACKS 

( 4  1 

=  14 

BOTH  SHIPS  SUNK  (]) 

=  7 

BOTH  SHIPS  SUNK  (2) 

=  4 

BOTH  SHIPS  SURVIVE 

(  3) 

=  0 

BOTH  SHIPS  SURVIVE 

14) 

=  0 

BOTH  SHIPS  SURVIVE 

(5) 

=  2 

OWN  SHIP  SUNK  (1) 

=  0 

CASE  II 

SUCCESSFUL  ATTACKS 

(5) 

=  33 

SUCCESSFUL  ATTACKS 

(b) 

=  180 

SUCCESSFUL  ATTACKS 

(7) 

=  149 

SUCCESSFUL  AT IACKS 

IS) 

=  47 

SUCCESSFUL  ATTACKS 

(9; 

=  57 

BOTH  SHIPS  SUNK  (b> 

=  24 

BOTH  Sll  1  PS  SUR\  1  VE 

(.7  ) 

=  ]n: 

BOTH  SIMPS  SURVIVE, 

(S) 

=  91 

BOTH  SHIPS  SURVIVE 

<9) 

=  10 

BOTH  SHIPS  SURVIVI 

(10) 

=  lb 

PCDET  1 


4lJ  . 


BOTH  SHIPS  SI  RV  1  V  K  lllj 


11 


BOTH  SHIPS  SIA'K  (L2)  =15 
BOTH  SHIPS  SUNK  (13)  =31 
OWN  SHIP  SUNK  (2)  =109 
OWN  SHIP  SUNK  (3)  =4 


TOTAL  SUCCESSFUL  ATTACKS 

TOTAL  BOTH  SHIPS  SUNK 

TOTAL  BOTH  SHIPS  SURVIVE 

TOTAL  OWN  SHIP  LOSS,  TARGET  SURVIVES 

TOTAL  OWN  SHIP  DETECTED  BEFORE  ATTACK 

EXCHANGE  RATIO  =  3.35  TO  1 


=  569 
=  81 
=  237 
=  113 
=  884 


30. 


APPENDIX  E 


COMPUTER  PROGRAM  FOR  THE  LINEAR  STATISTICAL  MODEL 


51. 


N1  *4 
N2  =  9 
SUM  X  =  0 

on  10  1=1. N2 

REA015.11)  Y(I)»(X(I»J),J=1,N1) 

11  FORMAT l 10F4 . 2 ) 

10  CONTINUE 

DO  20  1*1, Nl 
DO  21  J= 1  *  N2 
XT ( I  ,  J)  =  X( J, I ) 

21  CONTINUE 
2C  CONTINUE 
OC  1  1  =  1  *  N 1 
00  3  11=1, Nl 
00  2  J=1,N2 

SUMX=SlJMX+(  XT(  I  ,  J)  *X<  J,I  1)  ) 

2  CONTINUE 

xt x ( i  ,m*suMx 
SUMX=0 

3  CONTINUE 
1  CONTINUE 

CALL  GAUSS3C  N.EPS.XTX.XTXI ,KFR,9) 
IF(KER.NE.2)G0  TO  61 
WRITE<6,62) 

62  F0RMATU5X,8HSINGULAR) 

STOP 

61  SUMX=0 

00  50  J=1,N2 

WRITE(6.70)  (XT( I,J),I=l,Nl) 

70  F0RMAT(5X,4F8.4//I 
50  CONTINUE 

00  51  1*1, Nl 

WRITE! 6.71 )  (XTX(I.J)  ,J  =  l,Nl) 

71  F0RMAT(5X,4F13.3//) 

51  CONTINUF 

00  52  I  =  1 . N l 

WRITE<6,72)  (XTXIII,J),J=1,N1) 

72  F0RMAT(20X,4F20. 10//) 

52  CONTINUE 

00  4  12=1, Nl 
DO  5  13=1, N2 
00  6  14=1, Nl 

SUMX=SUMX*< XTXI ( I2,I4)*XTC 14,13) ) 

6  CONTINUF 

STO( I  2 , 1  3 )  =  SUMX 
SUMX=0 
5  CONTINUE 

4  CONTINUE 
SUNX=0 

DO  7  15=1, Nl 
DO  8  16=1, N2 

SUNX=SUNXM  STO<  15, T6)*YII6) ) 

8  CONTINUE 
8H ( I  5 ) =SUNX 
SUNX=0 

7  CONTINUE 
YTY=0 
BXY  =  0 


. 


52. 


00  30  I = 1 , M2 
yty=yty*i y(  imini 
3C  CONTINUE 
SUMX=0 

00  31  J= 1 , N2 
00  32  1=1, N1 

SUMX=SUMX+!8H< I ) ♦XT ( I ,J)I 

32  CONTINUE 
XSt J)=SUMX 
SUM  X=0 

31  CONTINUE 

00  33  1=1, N2 
BXY=BXY*(  XS<  I )*Y!II ) 

33  CONTINUE 
SUMX=YTY-BXY 
0=N2-N1 

SH  =  < SUMX ) /( 0  ) 

00  34  1=1, N1 

SIG(  I  )  =  (  XTXI  (  I,m*(SH) 

34  CONTINUE 
WRITEI6, 12) 

12  FORMAT!25X,8HBFTA  HA T , 2 5X , 9HS I GMA  HAT//) 

00  25  J= 1 , N 1 

MR  I T  E !  6  , 26) 8H  ( J) ,SIG( J) 

26  FORMAT! 15X,F30. 15, 5X,F30. 15// ) 

25  CONTINUE 

00  42  1=1, N2 
SUMX=0 

DO  41  J= 1 , N 1 

SUMX  =  SUM X+ ( RM( J ) *X ! I • J  > ) 

41  CONTINUE 
YHATU  )  =  SUMX 

YOIFC  I  )  =  Y(  I  l-YHATII) 

42  CONTINUE 

HR  I TE I  6 . 43 ) 

43  FORMAT! I 5X,4HY< I ) ,8X,7HYHATl I ) ,5X, 7HY0IF!  I  ) ) 
00  45  1=1, N2 

MR  I TE ( 6 • 44 )  Y(I)»YHAT(I)»YDIF!  I) 

44  FORMAT! I 3X,F  12.10, 2X , F 1 2 . 1 0 , 2X , F l ? . 1 0/ / ) 

45  CONTINUE 

00  80  1=1, N2 

REAO(5»81)Y! I) »(X<  t,J),J=l,Nl) 

81  FORMAT i 5  F4 . 2 ) 

8C  CONTINUE 

00  83  1=1, N2 
SUMX=0 

DO  84  J=  1  ,  N 1 
SUMX=SUMX+(RH< J)*X! I,J)  ) 

84  CONTINUE 

YH AT ( I )  =  SUMX 

YD  IF ( I )  =  Y! I ) -YH AT {  I  ) 

83  CONTINUE 
MRITE16.85) 

85  FORMAT! 15X,4HY( 1) , 8X »  7HYHAT ( I ) ,5X, 7HY0IF!  I  ) ) 
00  86  1=1, N2 

WRITE!6,87)  Y!  I  )  ,YHAT ( I ) , YOIF ( I) 

8  7  FORMAT! 13X, FI  2. 10, 2X , FI  2 , 10 , 2X , F 1 2 . 10/ /  ) 

86  CONTINUE 
MRITF<6,27) 

27  FORMAT  i 1H1 ) 

IOC  CONTINUE 

RETURN 

ENO 
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